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We present a systematic method for designing a perfect antireflection coating (ARC) for a semi-infinite 
one- dimensional (ID) photonic crystal (PC) with an arbitrary unit cell. We use Bloch wave expansion and 
time reversal symmetry, which leads exactly to analytic formulas of structural parameters for the ARC 
and renormalized Fresnel coefficients of the PC. Surface immittance (admittance and impedance) matching 
plays an essential role in designing the ARCs of ID PCs, which is shown together with a practical example. 
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Photonic crystals (PCs) have unique energy disper- 
sion due to coupling betweeB,.-»eKiodic materials and 
electromagnetic (EM) waves.liH'EHa In particular, the 
strong energy dispersion of the propagation modes of 
PCs has attracted much interest because it enables ap- 
plying PCs to add/drop multiplexers, dispersion co: 
pensators, polarization filters, and image processorsB 
These transmission-type applications require a. negligi- 
bly small reflection loss at the PC interfaced There- 
fore, applying effective antireflection interface structures, 
or antireflection coatings (ARCs) to the j»p«L-and out- 
put interfaces of the PCs is important .BH'EHO Several 
ARCs designed for the PCs have been reported on. The 
structural ideas of these ARCs weie-Jaased on the con- 
cepts of adiabatic intercapHectiora'LrEj or wave vector 
matching (k matching)t2Hlil for two-dimensional PCs. 
The reflectance at the interface of one-dimensional (ID) 
PCs has so far been calculated by plane-wave expan- 
sion, and, multiplication of the transfer matrix for a unit 
cell.E3lij'oE3lL3 However, these approaches do not tell 
directly the optimal ARC parameters, so the numerical 
calculation needs to be iterated until these parameters 
are optimized. In this letter, we derive analytic formulas 
of the structural parameters of an ARC that is applied 
to a ID PC. We deal with ID systems for simplicity, but 
the derived formulas would be applicable to multidimen- 
sional PCs under appropriate approximation. 

The structural parameters of a conventional ARC 
placed between two homogeneous media can be calcu- 
lated easily, if-. the refractive indices of the two media are 
known.E3E*ril In Fig. |l| (a), the ARC (region 2) with re- 
fractive index n 2 and thickness d is placed between two 
semi-infinite homogeneous media with refractive indices 
ni (region I) and 113 (region 3). These three regions are 
divided by two boundaries at z = and —d, where the z 
axis is defined as perpendicular to the surface of the ARC. 
Each region consists of linear and lossless dielectrics. The 
reflection coefficient at z = — d is then given by 
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FIG. 1. Antireflection coating (region 2) with thickness d 
and refractive index ri2 for (a) a semi-infinite homogeneous 
medium with refractive index nz and for (b) a semi-infinite 
ID PC with lattice constant A and periodic refractive index 
function n(z)(= n(z + A)). 



where k 2 is the normal component of the wave vector 
in the ARC, and 7*jj is the reflectkm, coefficient of the 
light propagating from region i to j.ExEj The reflection 
coefficient in Eq. ([!]) equals zero for normal incidence 
light when d = Ao/4n2|-and n<i = y^unjj, where Ao is the 
wavelength in vacuumO 

As in Fig. [l] (b), we replace the semi- infinite homoge- 
neous medium (region 3) in Fig. |l| (a) by a semi-infinite 
ID PC with lattice constant A and periodic refractive 
index function n{z). The function n(z) has an arbitrary 
spatial modulation in the unit cell. To determine the re- 
flection coefficient for this case, we expand the EM waves 
for the seaai-infinite PC by all the eigenmodes of the infi- 
nite PC.E3 These eigenmodes, which are extended into a 
complex k space (i.e., including the decaying waves), in 
the infinite ID PC paH-beiiObtained as the eigenvectors 
of a transfer matrixliJ'EJ'liata with a periodic boundary 
condition for a unit cell under a given frequency u>, a 
parallel component of a wave vector k|| , and a polariza- 
tion of light a. The reflection coefficient r of the electric 
field at the interface (z = —d) is determined by the con- 
tinuity of the tangential components of EM waves at two 
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boundaries (z = 0,—d). The expression of the derived 
reflection coefficient is formally the same as Eq. (]l]), ex- 
cept for f2,3. The reflection coefficient r 2 ,s for normal 
incidence light is modified into 



^2,3 (w) 



n 2 - N(w) 
n 2 +N(uj) 



(2) 



where N(oj) equals (u>) / e c, which is a normal- 

ized surface admittance of the Bloch waves. Note that 
in ±y£%(w), the "+"("-") sign applies to (Y^ 
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The surface admittances k f° r ^ ne ^ w0 orthogonal po- 
larizations are defined as 
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H k jE Kw , r k ( k } = H Kx /E kjy , at z = +0, (3) 



where k is the Bloch wave vector, and £?k,£(-ffk,$) with 
£ = x or j/ stands for the tangential component of the 
electric(magnetic) field of the propagating Bloch waves 
with a positive group velocity or the decaying Bloch 
waves in the positive z direction. These Bloch waves 
can be calculated from the forementioned transfer ma- 
trix, hence N(w) in Eq. (E) is obtained directly. In the 
derivation of Eqs. iO) and (|) we used the fact that time 
reversal symmetrycll inhibits the simultaneous appear- 
ance of the propagating and the decaying modes in the 
same direction for a given set {w, k|| and a} in a semi- 
infinite ID PC. Note that the function N(cJ) in Eq. (§) 
is generally complex due to the phase difference between 
-Ek;,£ and -Hk.r? (£> V = x or y) of the Bloch waves at the 
interface, however the imaginary part of N(u>) is zero if 
the surface of the semi-infinite PCs is a mirror plane in an 
infinite form of the PCs. Note also that multiple reflec- 
tions in the semi-infinite ID PCs, which are represented 
as plane waves, are renormalized into a single Fresnel 
coefficient in Eq. (g) by using the Bloch wave expansion. 

The pair values of refractive index n 2 and thickness 
d of the ARC for semi-infinite fD PCs that eliminate 
reflectance are determined, as follows. The extremal con- 
dition of the reflect anceEa R(— \r\ 2 ) with respect to d, 
i.e., dR/dd = 0, is written as 
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m = 0, 1, ■ 



(4) 



We call D the "normalized thickness" . By substituting d 
derived from Eq. (^) into Eq. (|f|) , we obtain a refractive 
index n 2 that eliminates the reflectance: 



ri2 = y/niRe(N)^ 



l n x - |7V| 2 /Re(A) 
n~i - Re(N) 



(5) 



The necessary and sufficient condition for the perfect 
ARC for semi-infinite ID PCs is D > 0, n 2 > 1 in Eqs. 
(Q) and (||), and d 2 R/dd 2 > under a given m > 1, 
N(u), and u. 

Next, we investigate the physical conditions encapsu- 
lated in Eq. (0). The condition Re(A r ) > is necessary 




FIG. 2. (a) Illustration of sufficient conditions in Eq. (6). 
The shaded areas indicate the region satisfying the conditions 
with — f3 2 — —0.08 (dotted line) as an example, (b) Magnifi- 
cation of region < a < 1 and —0.3 < / < of Fig. 2(a). 
Function / with rn = 2 is also plotted. 



to obtain a positive real number for n 2 , so the explicit 
form of n 2 > 1 with Eq.(^) is classified by using positive 
a = Re(N)/m: 



n 2 > 1 <^ 
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1 < a < oo 
a = 1 
< a < 1 



(6) 



where /(a, rii) = (a — l)(a — 1/nf), and (3 = lm(N)/ni. 

These sufficient conditions in the rhs of Eq. (ph are 
illustrated in Fig. [2] (a). As n\ changes from 1 to 00, 
the curve for / as a function of a varies from a solid one 
to a broken one continuously. The shaded area indicates 
the region satisfying the conditions of Eq. (^). When 
1 < a < 00, / > — /3 2 always holds because / > 0. The 
case a = 1 leads to (3 = due to / = 0. A magnification 
of region < a < 1 in Fig. || (a) is shown in Fig. || (b), 
where a curve for function / with n\ = 2 is added. The 
thick solid line shows the values of / and a that satisfy 
the third condition in the rhs of Eq. (^) for m = 2 and 
a given — f3 2 . This curve connects two crossing points 
(filled circles) of / and — (3 2 through the minimum point 
(open circle) of /. These three points are used to rewrite 
the inequality / < —f3 2 for < a < 1 in Eq. (||), which 
leads to two conditions: (1) the minimum point of / is 
less than or equal to — ft 2 , and (2) a is located between 
the two crossing points. 

By applying the previous, we can obtain the final form 
of the three sufficient conditions for n 2 > 1 in Eq. (||): 



(A) 1 < a < 00 , 

(B) a = 1 and (3 = , 

(C) < a < 1 , < g and 
-y/9 < 1 + 1/n? - 2a < y/g 



(7) 

(8) 

(9) 



where g((3,m) = (l — 1/nf) 2 — A/3 2 . Note that the con- 
dition (C) includes n\ 7^ 1 implicitly. Conditions (A)-(C) 
depend on n% and N(u>) only, so the n 2 needed to achieve 
ARCs for semi-infinite ID PCs can be directly obtained 
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Refractive Index n 2 Im[k] Re[k] Reflectance 

[Units of 2jt/A] 

FIG. 3. (a) Pairs of {yii and d} needed to achieve ARC for 
a semi-infinite ID PC. (b) Dispersion relation of propagating 
and decaying modes, (c) Reflectance with (solid line) and 
without (broken line) ARC. 



signing method would to a large extent be applicable 
to multidimensional PCs that have a position-dependent 
surface-immittance across the surface. 
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by using Eqs. (0)-(||) and d|). Accordingly, thickness d 
is also determined by Eq. (^ and d 2 R/dd 2 > 0. 

To give a practical example of ARCs for semi-infinite 
ID PCs, in Fig. ||(a) we plot the refractive index and 
the thickness needed to form the ARC for a semi-infinite 
ID PC. We assume that m = 1, that the unit cell of the 
PC consists of Si (n — 3.5) and SiC"2 (n = 1.5) with the 
same thickness A/2, and that the Si layer is the bound- 
ary layer between the semi-infinite PC and the ARC. In 
Fig. ||(b), the dispersion relation of the infinite ID PC is 
plotted and the band index is indicated by B1-B5. When 
the frequency changes along a band, the corresponding 
pair of [ri2(uj) and d(u>)] forms a continuous trajectory [as 
shown in Fig. 0(a)]. Note that in Fig. |^(a), only condi- 
tion (A) is needed to calculate ni and d, and frequency 
regions near the band edges that require n% > 4.5 are 
not plotted. Note also that the deviation in the normal- 
ized thickness from 1 is due to the phase shift between the 
electric and the magnetic fields at the interface (z = +0). 

In Fig. ||(c) we show the reflectance of the ID PC with 
(solid line) and without (broken line) an ARC. The ARC 
was designed at a selected frequency (ujK/2hc = 0.7), 
which is indicated by an arrow, where N(= 3.384— iO. 499) 
fulfills condition (A). The required values of ri2(= 1.868) 
and D{— 1.071) for the ARC are calculated from Eqs. (||) 
and (^) with m — 1, which is shown by a circle in Fig. 
0(a). In Fig. |](c), the ARC eliminates the reflectance 
(solid line) at the selected frequency (arrow). 

In conclusion, we have presented a systematic method 
for designing a perfect ARC for a semi-infinite ID PC. 
We derived exact formulas of structural parameters of 
the ARC and jEijesnel coefficient of the PC. The concept 
of k matchingE90 was not used for our derivations, and 
the Bloch wave vector in the infinite PC is not included 
explicitly in Eqs. (0) and (|^). The surface immittance 
(admittance and impedance) of the PC plays an essential 
role in the formulas. The adiabatic interconnectionljij'tij 
can be interpreted to be gradual immittance-matching 
structures. Under appropriate approximation, the de- 
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